Modeling micro-macro pedestrian counterflow in 
heterogeneous domains 

Joep Evers*t Adrian Muntean* 
November 8, 2010 

Abstract 

We present a micro-macro strategy able to describe the dynamics of 
crowds in heterogeneous media. Herein we focus on the example of pedes- 
trian counterflow. The main working tools include the use of mass and 
porosity measures together with their transport as well as suitable ap- 
plication of a version of Radon-Nikodym Theorem formulated for finite 
measures. Finally, we illustrate numerically our microscopic model and 
emphasize the effects produced by an implicitly defined social velocity. 
Keywords: Crowd dynamics; mass measures; porosity measure; social 
networks 
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1 Introduction 

One of the most annoying examples of collective behavioiQ is pedestrian jams 
- people get clogged up together and cannot reach within the desired time the 
target destination. Such jams are the immediate consequence of the simple ex- 
clusion process (181 124] , which basically says that two individuals cannot occupy 
the same position x E ft C R d at the same time t £ S :=]0, T[, where T g]0, oo[ 
is the final moment at which we are still observing our social network. 

Observational data (cf. e.g. [19]) clearly indicates that such jams typically 
take place in certain neighborhoods of bottleneck^] (narrow corridors, exits, 
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Figure 1: Schematic representation of the heterogeneous medium fl. The little 
black discs represent the pedestrians, while the dark gray zones are the parts 
where the pedestrians cannot penetrate (i.e. subsets of Q s ). The pedestrians are 
considered here to be the microscopic entities, while the grayish shadow indicates 
a macroscopic crowd; see Section |2 . 1 1 for the precise distinction between micro 
and macro made in terms of supports of micro and macro measures. 



corners, inner obstacles/pillars, ...). The effect of heterogeneities]^] on the overall 
dynamics of the crowd is what motivates our work. 

In this paper we start off with the assumption that inside a given room (e.g. 
a shopping mall), which we denote by i7, there are a priori known zones with 
restricted access for pedestrians (e.g. closed rooms, prohibited access areas, 
inner concrete structures Q whose union we call fl s . Let us also assume that 
the remaining region, say which is defined by f2 p := f2 — Q s , is connected. 
Consequently, f2 p is accessible to pedestrians. The exits of f2 - target that each 
pedestrian wants to reach - are assumed to belong to the boundary of fl p . The 
way we imagine the heterogeneity of is sketched in Figure [T] 

In this framework, we choose for the following working plan: Firstly, we 
extend the multiscale approach developed by Piccoli et al. [TU] (see also the 
context described in [5T] and [52]) to the case of counternowj^of pedestrians; then 
we allow the pedestrian dynamics to take place in the heterogeneous domain f2, 
and finally, we include an implicit velocity law for the pedestrians motion. The 
main reason why we choose the counterflow scenario [also called bidirectional 

3 Note that, for instance, Campanella et al. |8] give a different meaning to heterogeneity: 
they mainly refer to lack of homogeneity in the speed distributions of pedestrians. In [7] the 
geometric heterogeneities - obstacles - are introduced in the microscopic model. 

4 Note that some neighborhoods of these places can host, with a rather high probability, 
congestions! 

5 Two groups of people are moving in opposite directions. 
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flow [7J] out of the many other well-studied crowd dynamics scenarios is at least 
threefold: 

(i) Pedestrians counterflow is often encountered in the everyday life: at pedes- 
trian traffic lights, or just observe next week-end, when you go shopping, 
the dynamics of people coming against your walking direction [especially 
if you are positioned inside narrow corridors] . 

(ii) The walkers trying to move faster by avoiding local interactions with 
the oncoming pedestrians facilitate the occurrence of a well-known self- 
organized macroscopic pattern - lane formation; see, for instance [15j . 

(iii) We expect the solution to microscopic models posed in narrow corridors 
to be computationally cheap. Consequently, extensive sensitivity analy- 
ses can be performed and the corresponding simulation results can be in 
principle tested against existing experimental observations [TH1 H] ■ 

The presence of heterogeneities is quite natural. Pedestrians typically follow 
existing streets, walking paths, they trust building maps, etc. They take into 
account the local environment of the place where they are located. If the number 
of pedestrians is relatively high compared to the available walking space, then 
the crowd-structure interaction becomes of vital importance; see e.g. [5] for 
preliminary results in this direction. 

As long term plan, we wish to understand what are the microscopic mech- 
anisms behind the formation of lanes in heterogeneous environments. In other 
words, we aim at identifying links between social force-type microscopic models 
(see [UJ|5D], e -g-) an d macroscopic models for lanes (see [TH [3], e.g.) in the 
presence of heterogeneities. Here we follow a measure-theoretical approach to 
describe the dynamics of crowdfj^] Our working strategy is very much inspired 
by the works by M. Bohm [4] and Piccoli et al. [10] . 

The paper is organized as follows: In Section [2] we introduce basic model- 
ing concepts defining the mass and porosity measures needed here, as well as 
a coupled system of transport equations for measures. In Section [3] we present 
our concept of social velocity. Section [4] contains the main result of our paper 
- the weak formulation of a micro-macro system for pedestrians moving in het- 
erogeneous domains. We close the paper with a numerical illustration of our 
microscopic model (Section[5| exhibiting effects induced by an implicitly defined 
velocity. 



2 Modeling with mass measures. The porosity 
measure 

For basic concepts of measure theory and their interplay with modeling in ma- 
terials and life science, we refer the reader, for instance, to [T2] and respectively 

to [1 HQ no. 

6 The pedestrians are not exposed here to panic situations. 
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2.1 Mass measure 



Let fl C R d be a domain (read: object, body) with mass. Since we have in 
mind physically relevant situations only, we consider d £ {1,2,3}. However, 
most of the considerations reported here do not depend on the choice of the 
space dimension d. Let fi m (Q') be defined as the mass in CI' C O. Note that 
whenever we write SI' C Q, we actually mean that ST is such that SI' € B(fl), 
where B(Q) the a-algebra of the Borel subsets of ft. As a rule, we assume \x m 
to be defined on all elements of B(O). 

In Sections |2.1.1| and |2.1.2| we consider two specific interpretations of this 
mass measure that we need to describe the behavior of pedestrians at two sep- 
arated spatial scales. 



2.1.1 Microscopic mass measure 

Suppose that CI contains a collection of N point masses (each of them of mass 
scaled to 1), and denote their positions by {pk}k=i c ^> f° r N £N. We want 
fj, m to be a counting measure (see Sect. 1.2.4 in [T], e.g.) with respect to these 
point masses, i.e. for all SY € B(Cl): 

A*m(n') = #fe e n'}. (1) 

This can be achieved by representing p m as the sum of Dirac measures, with 
their singularities located at the pk, k £ {1,2,..., N}, namely: 

N 

Ihn = S P« ■ ( 2 ) 
fc=l 

We refer to the measure /x m defined by pi) as microscopic mass measure. 



2.1.2 Macroscopic mass measure 

Let us now consider another example of mass measure /x m . To do this, we 
assume that the following postulate applies to fi m : 

Postulate 2.1 (Assumptions on fj, m ). (i) fx m ^ 0. 



(ii) \i m is a-additive. 



(Hi) fj, m -C A d , where X d is the Lebesgue-measure in 



By Postulate 2.1 (i) and (ii), we have that fi m is a positive measure on Q, 
whereas (iii) implies that there is no mass pres ent i n a set that has no volume 
(w.r.t. \ d ). A mass measure satisfying Postulate ! 



2.1 



is in this context referred to 



as a macroscopic mass measure. Radon-Nikodym Theorerrj^] (see [T2] for more 

7 See [5] for a variant of this Theorem formulated for finite measures which is applied here. 
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details on this subject) guarantees the existence of a real, non-negative density 
p € L^ d (Sl) such that: 



M m (0') = / p(x)d\ d {x) for all ST £ B{tt). (3) 



Q' 



Similarly, we introduce time-dependent mass measures /it, where the time 
slice t G S enters as a parameter. 

2.2 Porosity measure 

Let SI C K d be a heterogeneous domain composed of two distinct regions: free 
space for pedestrian motion and a matrix (obstacles) such that fl = Q s U Sl p 
(disjoint union), where fl s is the matrix (solid part) of SI and Sl p is the free 
space (pores) . This notation is very much inspired by the modeling of transport 
and chemical reactions in porous media; see [2], e.g. 
Let /z p (f2') be the volume of pores in SI' C O. 

Postulate 2.2 (Assumptions on /i p ). (i) Hp ^ 0. 

(ii) /ip is a -additive. 
(Hi) Hp < A d . 



By Postulate 2.2 (i) and (ii), we have that Hp is a measure on f2. We refer 
to Hp as a porosity measure (cf. [1]). The absolute continuity statement in 
(iii) formulates mathematically that there cannot be a non-zero volume of pores 
included in a set that has zero volume (w.r.t. X d ). Assume that fl is such that 
A d (Sl) < oo. Then the Radon-Nikodym Theorem ensures the existence of a 
function <fi £ L+(ii) such that: 

H„(tf) = ( (j>d\ d for all O! £ B(fi). (4) 
Jo, 1 

Note that Hp(^') measures the volume of a subset of il' (namely of SI' n S7 p ). 
So, we get that 

H p (n') = \ d (n' n Sip) A d (Sl') for all Q' e B(n). (5) 

We thus have J n , (t)d\ d ^ f n ,d\ d , or / n ,(l - 0)dA d ^ 0. Since the latter 
inequality holds for any choice of Q' , it follows that 0^1 almost everywhere in 
SI. 

2.3 Transport of a measure 

For the sequel, we wish to restrict the presentation to the case d — 2. For 
our time interval S and for each i £ {1,2}, we denote the velocity field of the 
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corresponding measure by v l {t, x) with (t, x) G 5x0. Let also n\ , and juf be two 
time-dependent mass measures. Note that for each choice of i, the dependence 
on t of v l is comprised in the functional dependence of v l on both measures n\ , 
and /j 2 . This is clearly indicated in (|9j). The fact that here we deal with two 
mass measures fi] and fx 2 , transported with corresponding velocities v 1 and v 2 , 
translates into: 

'M + v .foV) = o, 

< d \ for all (t, x)eSxn. (6) 

t^ + v. ( „y) = o, 

These equations are accompanied by the following set of initial conditions: 

l4 = f4> ast = 0forie{l,2}. (7) 

It is worth noting that (|6| is the measure-theoretical counterpart of the Reynolds 
Theorem in continuum mechanics. To be able to interpret what a partial differ- 
ential equation in terms of measures means, we give a weak formulation of ([6]). 
Essentially, for all test functions ip 1 ,^ 2 € Cq (O) and for almost every t € S, the 
following identity holds: 

— / ip l (x)dnl(x) = [ v i {t,x)-Vtl) i {x)d^ t {x)iov alH e {1,2}. (8) 

Definition 2.1 (Weak solution of ([6])). The pair }t^o, {/^t }t^o) * s called a 
weak solution of |6j), if for a/He {1.2} £/ie following properties hold: 

1. the mappings t i— > J* n tp" 1 (x)d^x\{x) are absolutely continuous for all ip l £ 

2. ^eL 2 (S;tfj(fi)); 

5. Equation Q) is fulfilled. 

We refer the reader to [9] for an example where the existence of weak solu- 
tions to a similar (but easier) transport equation for measures has been rigor- 
ously shown. 

3 Social velocities 

We follow very much the philosophy developed by Helbing, Vicsek and coauthors 
(see, e.g. [15] and references cited therein) which defends the idea that the 
pedestrian's motion is driven by a social force. Is worth noting that similar 
thoughts were given in this direction (motion of social masses/networks) much 
earlier, for instance, by Spiru Haret [T3] and Antonio Portuondo y Barcelo [25] . 
Moreover, other authors (for instance, Hoogendoorn and Bovy [T5]) prefer to 
account also for the Zipfian principle of least effort for the human behavior. We 
do not attempt to capture the least effort principle in this study. 
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3.1 Specification of the velocity fields v l 

Until now, we have not explicitly defined the velocity fields v l (i g {1,2}). Very 
much inspired by the social force model by Dirk Hclbing et al. |14j , the velocity 
of a pedestrian is modeled as a desired velocity v des perturbed by a component 
v\ i 2i- The latter component is due to the presence of other individuals, both 
from the pedestrian's own subpopulation and from the other subpopulation. 
The desired velocity is independent of the measures /i* and fj, 2 , and represents 
the velocity that an agent would have had in absence of other pedestrians. 

For each i £ {1, 2}, the velocity if is defined by superposing the two velocities 
t)j oc and follows: 

v t (t,x):=v t dcs {x)+v t M ii 2 ] (x), for all t £ (0, T) and x £ VL. (9) 

For a counterflow scenario, the desired velocities of the two subpopulations 
follow opposite directions. We thus take 

«d»(aO = «dc S e R 2 

fixed (for i £ {1,2}) and 

v dcs = ~ u dcs- 

The component u^i ^ models the effect of interactions with other pedestri- 
ans on the current velocitjj^] Since the interactions between members of the 
same subpopulation differ (in general) from the interactions between members 
of opposite subpopulations, we assume that v] x 2 i consists of two parts: 

y-x 



+ r v {\y-x\)g{tf xy )^ M(y), (10) 

Jq\{x} \y-x\ 



for i £ {1,2}, where j = 1 if i = 1 and vice versa. In (10) we have used the 
following: 

• / own and / opp are continuous functions from M + to M, describing the effect 
of the mutual distance between individuals on their interaction. Compare 
the concept of distance interactions defined in [25] . / own incorporates the 
influence by members of the same subpopulation, whereas / opp accounts 
for the interaction between members of opposite subpopulations. / own is 
a composition of two effects: on the one hand individuals are repelled, 
since they want to avoid collisions and congestion, on the other hand they 
are attracted to other group mates, in order not to get separated from 
the group. / opp only contains a repulsive part, since we assume that 
pedestrians do not want to stick to the other subpopulation. 



'The interactions we are pointing at are nonlocal. 
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• (x % xy denotes the angle between y — x and v l dcs (x): the angle under which 
x sees y if it were moving in the direction of v dcs (x). 

• g is a function from [— 7r, tt] to [0, 1] that encodes the fact that an individ- 
ual's vision is not equal in all directions. 

Regarding the specific choice of / own , / opp and g we are very much inspired 
by [13] and [TU], e.g. However we do not use exactly their way of modeling 
pedestrians' interaction forces. We list here the following forms for the functions 
/ own , / opp and g that match the given characterization: 

if s < /?°pp- 

/° PP (s) := I * V S 2 (R° pp )2j> . (11) 




if s > R° pp , 



if s s? i?° wn ; 



/ own (s) := { \s R™ n J\s " 11 ' (12) 

if a > R°r, 

g(a) := a + (1 - o) 1 + c ° s ( a \ for a G [-tt.tt]. (13) 

Here F opp and F own are fixed, positive constants. The constants R° pp , R° wn 
(radii of repulsion) and i?° wn (radius of attraction) are fixed and should be cho- 
sen such that < R° wn < R° WD and < R° pp . Furthermore, the restriction 
max{i?° wn , i?° pp } < L has to be fulfilled. The interaction / opp is designed such 
that individuals "feel" repulsion (i.e. / opp < 0) from another pedestrian if they 
are placed within a distance R° pp from one another. The corresponding state- 
ment holds for f own if individuals are within the distance R° wn . Additionally, 
an individual is attracted to a second individual if their mutual distance ranges 
between i?° wn and R° wn . 

The function g ensures that an individual experiences the strongest influence 
from someone straight ahead, since ^(0) = 1 for any u £ [0, 1]. The constant a is 
a tuning parameter called potential of anisotropism. It determines how strongly 
a pedestrian is focussed on what happens in front of him, and how large the 
influence is of people at his sides or behind him. 

In the remainder of this section, we suggest four different alternatives for the 
definition of ^ by indicating various special choices of distance interactions 



and visibility angles (conceptually similar to a l xy ) as they arise in (10 1. All 

of them boil down to including an implicit dependency of the actual velocity 

v l = wL, + ff i 2i- Note that this effect increases the degree of realism of the 

dcfa KXJ . .... 

model, but on the other hand it makes the mathematical justification of the 

corresponding models much harder to get. 
3.1.1 Modification of the angle a l xy 

We defined the angle a xy as the angle between the vector y — x and v des (x). 
However this is not a good definition if the pedestrian in position x is not moving 
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in the direction of v^ os (x) (or, in a broader sense, if the actual speed cannot be 
approximated sufficiently well by the desired velocity) . Therefore we suggest to 
define a xy = a xy (t) as the angle between y — x and v l (t, x). 

3.1.2 Prediction of mutual distance in (near) future 

Up to now the functions f own and / opp depended on the actual distance between 
x and y at time t. However pedestrians arc likely to anticipate on the distance 
they expect to have after a certain (small) time-step (say, some fixed At € R). 
In practice, this means that at a time t £ S a person will modify his velocity 
(either in direction, or in magnitude, or both) if he foresees a collision at time 
t + At e S. 

To predict the mutual distance between x and y at time t + At, the current 
velocities at x and y are used for extrapolation. The predicted distance is: 



\(y + v(y,t)At) — (x + v(x,t)At)\. Consequently, sticking to the notation in (10), 
the interaction potential / own and / opp should depend on \(y+v i (t, y)At) -(x + 
v t (t,x)At)\ and on \ (y + v j (t 7 y)At) — (x + v l (t,x)At) \ respectively (where j = 1 
if i = 2 and vice versa). 

3.1.3 Prediction of mutual distance within a time interval in the 
(near) future 

The disadvantage of using \ (y + v(y, t)At) — (x + v(x, t)At)\ is that At is fixed. 
A pedestrian can thus only predict the distance at an a priori specified point 
in time in the future. However, people are able to anticipate also if they expect 
a collision to occur at a time that is not equal to t + At. We assume now that 
we are given a fixed Ai max 6 M + such that an individual can predict mutual 
distances by extrapolation for any time r £ (i, i + Ai max ). Thus, Ai max imposes 
a bound on how far can an individual look ahead into the future. To capture 
this effect, we suggest to replace f° wn (\y — x\) and / opp (|y — x\) by: 



1 



At r , 



and 



At mayL 



f° wn (\(y + v l (t, y)r) -(x + v%x)T)\)dr, (14) 
/ opp (| (y + v j (t, y)r) -(x + v*(t, x)r)\)dT, (15) 



respectively. 



3.1.4 Weighted prediction 

Since an individual probably attaches more value to his predictions for points 
in time that are nearer by than others, one additional modification comes to 
our mind. Let h : [t, t + Ai max ] — > [0, 1] be a weight function. Then instead of 



(14) and (15), we propose 

At max 

r™(|(y + v*(t, y)r) - (x + v\t, x)r)\)h{r)dr, (16) 



1 
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and 

f ™ X f°PP(\(y + V 3( t ,y)T)-(x + V l (t,x)T)\)h(T)dT. (17) 



max JO 



If h is decreasing, then the influence of t\ is larger than the influence of ti, if 
t\ < t-2 (which matches our intuition). 

3.2 Two-scale measures 

We now consider the explicit decomposition of the measures p\ and /j^. Let the 
pair (0i, 9 2 ) be in [0, l] 2 , and consider the following decomposition of y\ and /x 2 : 

f4 = e i r4 + (l-9i)Mi, i€{l,2}. (18) 

Here, mj is a microscopic measure. We consider {p\.(t)}k=i C £1 to be the 
positions at time t of N l chosen pedestrians, that are members of subpopulation 
i. We want m\ to be a counting measure with respect to these pedestrians, i.e. 
for all Q' <= B(fi): 

mi(n') = #{p£(t) € fi'}, ie{l,2}. (19) 



We thus define ml as the sum of Dirac masses (cf. Section 2.1.11, centered at 
thep* fe , k= l,2,...,iV 4 : 



TV* 

m * = EV fc (t)- *e{i,2}. (20) 
k=l 

Ml is the macroscopic part of the measure, which takes into account the part 
of the crowd that is considered continuous. We consequently have Ml <C A 2 , 
since a set of zero volume cannot contain any mass. Note that we are thus 
in the setting of Section |2.1.2| Now, Radon-Nikodym Theorem guarantees the 
existence of a real, non-negative density p l (t, •) 6 L^ 2 (^) such that: 



M l t (n') = / p\t,x)d\ 2 (x) (21) 

JSl' 

for all Q! e B{Q) and all i € {1,2}. 

4 Micro- macro modeling of pedestrians motion 
in heterogeneous domains 

We have already made clear that we want to model the heterogeneity of the 
interior of the corridor. In practice this means that pedestrians cannot enter 



all parts of the domain. As described in Section 2.2 we have a measure n P 



corresponding to the porosity of the domain (which is fixed in time) . However, 



we note that the concept of porosity (cf. Section 2.2) is a macroscopic one 
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For this reason only the macroscopic part of the mass measure in ( 18 ) needs 
some modification with respect to the porosity. In this context, one should 
be aware of the analogy with mathematical homogenization. This technique 
distinguishes between microscopic and macroscopic scales, where we also see 
that some (averaged) characteristics are only defined on the macroscopic scale. 
For more details, the reader is referred to [5] or [T7]. In K 2 , we have fi p « ) 2 . 
Furthermore Ml <C p p for i G {1, 2} and a.e. t £ S. This is obvious, since no 
pedestrians can be present in a set that has no pore space (i.e. zero porosity 
measure). A basic property of Radon-Nikodym derivatives now gives us: 

dMl dMidu v r „ „ , 

— i = -—^-^ i G {1, 2} for almost every t € S. (22 
dX 2 dp p dX 2 

dAl i dp 

We have already defined p l (t,-) :— — and <j) :— — |. If we now denote 

d\ 2 d\ 2 

dM l 

by p l (t,-) the Radon-Nikodym derivative the following relation holds: 

d/i p 

p l {t,-) = p l (t,-)4>(-) for all i £ {1,2}. 

4.1 Weak formulation for micro-macro mass measures 

We now have the following measure: 

l4 = e t m\ + (1 - 6i)Mi, i G {1, 2}, (23) 



as was given in (18), where now: 



H — 

k=l 



^S pm , dMi(x)=p i (t,x)<P(x)d\ 2 (x). (24) 



This specific form of the measure will now be included in the weak formulation 
([8|, with velocity field ([£])-( 10). The real positive numbers 0j (i G {1,2}) are 



intrinsic scaling parameters depending on N l 

The transport equation (Tsl) takes the following form: 



, N' . 
^(^^^(PU*)) + (!-»<) / ^(x)p%x)<f>(x)d\ 2 (x) 
k=l n 

0i y2v*(t,pl(t)) ■ VtfipUt)) + (1 - 0i) / v%x) ■ V^(x)p\t,x)^(x)dX 2 (x), 

(25) 

for all i G {1,2}. Here we have used the sifting property of the Dirac distribu- 
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tion. In the same manner, we specify v\ ul u?1 from ( 10 ) as 

p%(t)-x 



\pW)-A 

-(1-610 / /^(ly-xDsK^^^pX^y)^)^ 2 ^) 



fc=i 



+(1-0,-) / / opp (|y-*|)<7«)^ »Mtf)dV(iO, 

Jn W ~ x \ 

for i 6 {1, 2}, and j as before (j = 1 if i = 2 and vice versa). We have omitted 
the exclusion of {x} from the domain of integration (in the macroscopic part), 
since {a;} is a nullset and thus negligible w.r.t. A 2 . Note that the sums may be 
evaluated in any point x G fi (not necessarily x — p\.{t) for some i and k)\ the 
integral parts may also be evaluated in all x, including x = p l k (t) for some i and 
k. 



5 Numerical illustration 

We wish to illustrate now the microscale description of a counterflow scenario 
(i.e. for 0i = 6 2 = 1) by presenting plots of the configuration of all individuals 
situated in a given corridor at specific moments in time. 

We consider a specific instance in which there are in total 40 individuals (20 
in each subpopulation) . The dimensions of the corridor are d — 4 and L = 20. 



The velocity is taken as defined in ( 10 1-( 13 ) . Furthermore, the following model 



parameters are used: v\ cs = 1.34e 1; ^ cs = -1.34e 1; F opp = 0.3, F own = 0.3, 
i?° pp = 2, R° wn = 2, R° wn = 0.5, F w = 0.5, R w = 0.5, a = 0.5. 

In Figure[2] we show the configuration in the corridor at times t = 0, t = 7.5, 
and t — 15. The individuals of the subpopulation 1 are colored blue, while the 
individuals of the subpopulation 2 are colored red. Clearly, self-organization 
can be observed in the system: Pedestrians that desire to move in the same 
direction form lanes (in this case, three of them). This feature is observed and 
described extensively in literature, cf. e.g. [T5] . 

Another feature, pointed out by Figure [2j is the following: Within the three 
already formed lanes, small clusters of people are formed. This flocking is a 
result of the typical choice for / own in ( |l~2"| . Members of the same subpopulation 
are repelled if their mutual distance is in the range (0, i?° wn ); they are attracted 
if their mutual distance is in the range (R° wn , R° wn ). No interaction takes place 
if individuals are more than a distance R° wn apart. 

The attraction part of the interaction causes individuals that are already rel- 
atively close to get even closer, until they are at a distance R° wn . For distances 
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Figure 2: The simulation of a crowd's motion in a corridor of length L = 20 
and width d = 4. Each of the two sub-populations consists of 20 individuals. 
The images were taken at t — (left), t = 7.5 (middle), t = 15 (right). 

around R° wn , there is an interplay between repulsion and attraction, eventu- 
ally leading to some equilibrium in the mutual distances between neighboring 
individuals in one cluster. In Figure [2j we observe self-organized patterns even 
within clusters. 
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